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DYNAMICS  OF  EVAPORATION  OF  Csl  MOLECULES  FROM 
INDIVIDUAL  [Cs(CsI)„]^  CLUSTERS 

HYUN  JIN  HWANG^  and  M.  A.  EUSAYED 
Department  of  Chemistry  and  Biochemistry, 

University  of  California 
Los  Angeles,  California  90024,  USA 


Abstract 


A  brief  review  is  first  ^ven  of  the  different  types  of  clusters  and  their  synthesis.  The 
results  on  the  determination  of  the  mechanisms  of  the  evaporation  of  one  and  two  Csl 
molecules  from  individual  mass  selected  [Cs(CsI)„]'*^  clusters  (made  by  sputtering 
techniques)  are  discussed  in  terms  of  their  observed  stabilities,  their  observed  aiverage 
kinetic  energy  released  and  its  distribution  during  the  evaporation  process.  It  is  found 
that  the  stability  towards  evaporation,  rather  than  formation  rates,  is  the  ori^  of 
observing  clusters  with  mag^c  numbers.  Furthermore,  comparison  with  theoretical 
models  showed  that  the  evaporation  process  of  Csl  molecules  from  these  clusters  is 
statistical.  Finally,  it  is  shown  that  the  evaporation  of  two  Csl  molecules  mvolves  the 
fission  of  a  dimer  rather  than  a  sequential  evaporation  of  two  Csl  monomers. 
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I.  Introduction;  Cluster  Types,  Synthesis  and  Detection: 

It  is  usually  accepted  that  the  study  of  clusters  gives  us  a  systematic 
understanding  of  the  physical  and  chemical  properties  of  matter  as  a  function  of  its 
size.  It  is  hoped  that  questions  such  as  how  many  atoms  one  needs  to  have  a  metal 
or  how  many  molecules  are  required  to  have  the  properties  of  a  molecular  liquid  or 
solid  can  be  answered.  Of  course,  the  answer  will  depend  on  the  length  or  time  scale 
of  the  experiment  used  to  measure  such  a  property  as  a  function  of  size. 

Clusters  are  classified  according  to  the  type  of  bonding  believed  to  hold  the 
cluster  together  or  the  type  of  entity  that  leads  to  cluster  formation.  Van  der  Waals 
clusters^  are  formed  between  rare  gas  atoms  or  stable  neutral  molecules  or  between 
rare  gas  atoms  and  stable  molecules.  Small  van  der  Waals  clusters  (in  particular 
dimers)  have  received  a  great  amount  of  attention  in  the  field  of  molecular 
dynamics,^  as  they  represent  an  interaaion  between  these  entities  during  "half 
collisions.  The  flow  of  eneigy  between  intramolecular  vibrations  and  the  motion 
along  the  van  der  Waals  dissociation  coordinate  has  been,  and  is,  of  great  theoretical 
and  experimental  interest.  They  have  recently  been  used  to  define  zero  time  in  rapid 
biomolecular  photochemical  reactions^  between  the  cluster  components  using 
femtosecond  spectroscopy.  Van  der  Waals  clusters  have  been  studied  by  mass 
spectrometry,  infrared,  Raman  and  optical  spectroscopy.  Clusters  in  which  hydrogen 
bond  might  play  a  role  in  the  binding,  e.g.,  those  involving  organic  acids,  or  sdcohols, 
NH3,  H2O,  are  termed  H-bonded  clusters."* 

Clusters  made  of  metal  atoms  are  termed  metallic  clusters.  Here  the  name  is 
based  on  the  fact  that  atoms  resulting  from  evaporated  metals  form  the  cluster. 
Certainly,  the  bonding  in  small  metallic  clusters  is  not  exactly  of  the  metallic  type 
found  in  bulk  metals.  For  alkali  and  noble  metals,  the  "shell"  or  Jellium  model  is 
found^  to  explain  the  observed  clusters  with  magic  numbers  (those  with  extra  stability 
compared  to  the  ones  following  them  in  size).  Recent  studies^  in  mercury  clusters 
(Hgjj),  suggest  that  for  n  s  13  van  der  Waals  forces  are  involved,  at  30  s  n  <  70 
covalent  bonds  are  formed  and  for  n  ^  100  metallic  bonding  takes  over.  A  large 
number  of  studies  on  physical  and  chemical  reactivity  of  transition  metal  clusters  has 
been  carried  out.^’®  Detailed  structural  calculations  for  large  transition  metal  clusters 
have  not  yet  been  attempted.  Due  to  the  large  number  of  unpaired  electrons,  these 
clusters  possess  a  large  number  of  electronic  states  that  are  sufficiently  close  in 
energy  to  allow  too  many  structures  (isomers)  of  comparable  energy, 

Nonmetallic  clusters,^  those  made  of  nonmetals,  e.g.,  carbon,  silicon,  are 
covalently  bonded.  Active  work  was  carried  out  by  Smalley  and  Kroto^®  and  their 
groups  to  study  a  cluster  which  is  expected  to  have  a  structure  with  special 
stability  and  is  thus  observed  with  high  intensity  in  the  mass  spectrum  of  carbon  ion 
clusters.  A  solid  which  is  mostly  and  some  Ojq  has  recently  been  isolated^^  in 
the  laboratory  by  condensing  the  vapor  formed  from  the  sublimation  of  carbon  rods 
in  helium  atmosphere  by  passing  high  electric  currents  through  it.  Studies  of  clusters 
made  from  semiconducting  material,^  e.g.,  CaAs  have  also  been  carried  out. 

Clusters  made  from  ionic  solids,  e.g.,  alkali  halide  salts,  have  also  been 
studied, both  experimentally  and  theoretically.  Clusters  with  magic  numbers  are 
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found  to  fit  a  lattice  with  sides  having  integer  number  of  ions.  Ousters  of  equal 
number  of  cations  and  anions  (i.e.,  neutral  clusters)  or  with  excess  cations  (e.g., 
[M(MX)^]''‘),  with  excess  anions  [X(MX)j,]‘^  or  with  excess  electrons^  [M^X^.jJ.e 
have  all  been  observed  and  studied. 

Over  the  years,  a  number  of  techniques  have  been  used  to  synthesize  clusters. 
In  solution,  metallic-containing  clusters  are  synthesized  by  use  of  colloidal  and 
electrochemical  techniques.^^  Metal  clusters  are  formed  from  their  polynuclear 
complexes  by  different  ligand  stripping  methods.^  Metal  ion  clusters  have  also  been 
made  in  liquid  metal  ion  sources^”  by  applying  a  high  electric  field  to  a  tungsten  tip 
dipped  in  the  liquid  metal  of  which  cluster  ions  are  to  be  made. 

Clusters  are  made  by  aggregation  on  a  support,  e.g.  in  the  rare  gas  matrix 
isolation  method^^  or  on  different  types  of  support,  e.g.,  molecular  sieve  host  and 
functionalized  polymers.^*  Matrix  isolation  methods  are  shown  to  be  especially 
useful  in  the  study  of  the  structure  of  very  small  clusters  using  spectroscopic  methods. 

A  number  of  methods  have  been  used  to  synthesize  gaseous  clusters,  by 
continuous  effusive  expansion  of  oven-heated  vapors  of  elements  in  the  pure  phase 
or  in  rare  gases^^  (oven  gas  aggregation  techmque^).  Pulsed  supersonic  cooling 
expansion  techniques  in  rare  gases  (He  or  Ar)  have  been  most  useful  for  the  making 
of  many  of  the  weak  van  der  Waals  or  H-bonding  clusters.  The  pulsed  mode  allows 
the  use  of  high  backing  pressures  of  the  rare  gas  used  and  thus  leads  to  much  lower 
temperatures  for  the  clusters  formed.  The  mole  fraction  of  the  gas  or  vapor  of  which 
the  clusters  are  to  be  made  is  adjusted  so  as  to  allow  for  collisions  between  them  to 
take  place  as  they  cool  off  during  the  expansion  process.  Dimers,  trimers  and  higher 
clusters  as  well  as  clusters  with  the  rare  gas  atoms  used  are  formed  under 
appropriate  expansion  conditions.  For  nonvolatile  material,  pulsed  laser  ablation, 
with^^  or  without^^  supersonic  expansion  of  the  vapor  in  rare  gases,  is  used  to 
produce  clusters  and  their  positive  and  negative  ions.  The  use  of  supersonic 
expansion  in  conjunction  with  pulsed  laser  ablation  has  greatly  increased  the  size  and 
types  of  the  clusters  made. 

Another  technique  used  for  making  gaseous  clusters  from  nonvolatile  material 
is  by  use  of  the  sputtering  technique  in  which  high  energy  particles,  e.g.,  electrons, 
rare  gas  cations  or  atoms,  are  used  for  heating  the  material.  Next  to  the  oven 
technique,  making  clusters  from  RF  spark  sources^  followed  shortly  in  time  the  use 
of  ovens.  Rare  gas  ion  sputtering^"^  and  now  fast  atom  bombardment  (FAB)^^  using 
netural  xenon  are  very  useful  in  studies  of  alkali  halides  and  high  boiling  metals. 

The  most  common  ways  of  determining  the  cluster  distribution  and  its  changes 
upon  reaction  are  quadrupole^  and  time-of-flight  mass  spectrometers.^  If  a  pulsed 
laser  is  used  for  ionization,  time-of-flight  is  the  spectrometer  of  choice  since  the 
whole  mass  spectrum  is  recorded  for  each  laser  shot.  This  eliminates  the  laser  pulse- 
to-pulse  intensity  fluctuations.  Cluster  reactivity  and  its  size  dependence  is  best 
studied  by  use  of  ion-cyclotron  resonance  mass  spectrometry.^ 

In  this  article  we  make  a  detailed  study  of  the  evaporation  dynamics  of  one  and 
two  Csl  molecules  from  mass  selected  lCs(CsI)^]'‘^  ionic  cluster  ions.  From 
measurement  of  evaporation  probability,  kinetic  energy  released  and  its  distribution 
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upon  evaporation  as  a  function  of  cluster  size,  the  following  conclusions  are  reached: 
1)  the  rate  of  evaporation  from  the  clusters,  rather  than  the  rate  of  their  initial 
formation  from  sputtering  techniques,  is  the  origin  of  the  appearance  of  clusters  with 
magic  numbers  in  the  mass  spectrum,  2)  the  evaporation  process  is  statistical,  and 
3)  the  loss  of  two  Csl  molecules  is  found  to  result  from  a  fission  of  a  cesium  iodide 
dimer  rather  than  a  sequential  loss  of  two  Csl  molecules.  An  account  of  this 
observation  is  given. 

II.  DYNAMICS  OF  EVAPORATION  OF  Csl  MOLECULES  FROM  [CsCCsUn]"^ 
CLUSTERS:^’’^® 

Alkali  halide  cluster  ions  of  the  form  [M(MX)jj]‘'’  have  been  extensively  studied 
both  theoretically^^  and  experimentally.^'  T^e  ease  of  their  formation  by  fast  atom 
or  ion  bombardment^  (by  sputtering  methods)  have  made  them  attractive  to  study 
experimentally.  The  simple  polarizable  ionic  bonding  model  is  found  to  yield  good 
values  for  their  formation  energies.'^^ 

A  strong  size  dependence  was  observed  for  the  cluster  ^eld  observed  in  the 
mass  spectrum  of  the  clusters  made  from  sputtering  methods. For  example,  for 
clusters  of  the  composition  [M(MX)J'*^,  mass  peaks  at  n  =  6,  9,  13,  22,  31,  37  ... 
were  found  to  have  unusually  high  intensities  compared  to  those  following  them  in 
mass.  In  earlier  studies,^^^^  there  was  an  attempt  to  attribute  the  high  yield  to  larger 
rates  of  formation  of  these  clusters  during  the  sputtering  process.  However,  Ens  et 
al.^^  showed  that  the  mass  spectrum  obtained  shortly  after  the  cluster  formation  (as 
opposed  to  the  much  longer  time  scale  of  observation  in  the  conventional  method 
using  a  sector-field  mass  spectrometer)  has  no  magic  numbers,  i.e.,  the  mass  peak 
intensity  is  found  to  decrease  exponentially  with  size. 

In  our  studies,  we  have  moved  away  from  a  study  of  the  total  mass  spectrum. 
We  examine  the  evaporation  pattern  of  size-selected  [Cs(CsI)n]'*^  after  a  flight  of  a 
fixed  distance.  The  evaporation  probability  during  this  time  is  determined  as  a 
function  of  the  cluster  size.  Furthermore,  by  using  a  high  resolution  sector-field  mass 
spectrometer,  the  kinetic  energy  released  during  the  evaporation  process  as  well  as 
its  distribution  are  determined  for  the  evaporation  of  one  or  two  Csl  molecules.  A 
number  of  questions  we  have  attempted  to  answer  in  our  studies  are:  a)  Can 
variation  in  the  evaporation  probability  with  size  account  for  the  observed  clusters 
with  magic  numbers?  b)  Is  the  evaporation  process  statistical  and  independent  of 
size?  c)  Does  the  loss  of  two  Csl  molecules  occur  sequentially  or  as  a  result  of  a 
dimer  fission? 

A.  ExperimentalP'^ 

We  have  used  a  VG  analytical  ZAB-SE  mass  spectrometer  fitted  with  a  fast 
atom  bombardment  gun  (model  FAB  1  IB)  in  the  single  focus  mode.  Cluster  ions  are 
produced  with  7  keV  Xe  atom  bombardment  of  Csl  target  prepared  by  evaporating 
aqueous  salt  solution  on  Ni  foil.  Positively  charged  cluster  ions  are  accelerated  to 
8.1  keV  through  a  two  stage  acceleration  plate  and  analyzed  by  a  66  cm  radius  35” 
magnetic  sector.  The  observed  mass  spectrum  showed  sharp  mass  peaks 
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corresponding  to  the  parent  cluster  ions  Cs(CsI)°''’,  and  diffuse  mass  peaks 
corresponding  to  the  daughter  cluster  ions  formed  mainly  by  evaporation  of  one  or 
two  Csl  molecules  in  the  field  fi’ee  drift  region  (between  the  acceleration  region  and 
the  magnetic  sector).  Evaporation  pathways  of  the  (size-selected)  parent  clusters 
[Cs(CsI)“]'''  with  1  s  n  s  17  were  determined  fi'om  the  observed  nominal  mass  values 
m^  of  the  daughter  cluster  ions  [Cs(CsI)“]'*'  with  m=n-l  or  n-2,  by  using  the  relation 
m  =  where  m^  and  m2  are  the  actual  mass  values  of  the  parent  and 

daughter  cluster  ions.  To  measure  the  kinetic  energy  release  distribution  as  well  as 
the  evaporation  probability  for  each  evaporation  pathway,  spectra  of  the  daughter 
and  corresponding  parent  ions  were  taken  by  scanning  the  acceleration  voltage 
after  tuning  the  magnetic  sector  field  to  the  mass  value  of  each  ion  at  the 
acceleration  voltage  of  8.1  kV.  During  the  experiments,  the  pressure  is  kept  at 
~2xl0'*  Torr  in  the  acceleration  region  and  at  -5x10'®  Torr  in  the  analyzer  or  field- 
free  region.  Thus  no  collisions  occur  during  the  evaporation  process.  Under  our 
experimental  conditions,  one  calculates  that  the  cluster  with  n  =  13  spends  57.1  ps  in 
its  flight  through  the  field-free  region. 

B.  Clusters  with  Magic  Numbers:  Variation  in  Evaporation  Probability: 

The  mass  peaks  corresponding^^  to  n  =  13,  22,  37,  52,  62,  are  found  to  be 
followed  by  weaker  mass  peaks,  suggesting  special  stability  for  clusters  having  these 
values  of  n,  the  magic  numbers.  Interestingly  enough  the  numbers  of  atoms  in  these 
magic  number  clusters  agree  with  "cube-like"  structures  with  edges  of  k  x  f  x  m, 
where  k,  f,  and  m  are  integers.  This  suggests  a  rectangular  or  cubic  ionic  lattice  that 
maximizes  the  electrostatic  interaction  energy  per  atom  (ion)  for  the  clusters  with 
magic  numbers.  Of  the  clusters  with  magic  numbers,  those  having  cubic  lattices  with 
k=f=m  are  found  to  have  more  intense  mass  peaks  [e.g.,  those  for  n=13  (3x3x3) 
and  62  (5x5x5)].  All  these  observations  are  in  agreement  with  the  polarizable-ion 
model  calculations.^^ 

The  question  that  arises  is  whether  this  unusual  stability  manifests  itself  in  an 
increased  rate  of  formation  of  these  clusters  during  the  sputtering  process,  or  if  the 
probability  of  formation  decreases  smoothly  with  the  cluster  size,  thus  allowing  the 
variation  in  the  evaporation  probability  to  be  the  origin  of  clusters  with  magic 
numbers.  Indeed  Ens  et  al."  have  demonstrated  that,  initially,  the  cluster  yield 
decreased  smoothly  with  size.  Below  we  show  that  the  variation  in  the  probability 
of  evaporating  one  or  two  Csl  molecules  can  lead  to  mass  intensity  distributions 
showing  the  magic  number  characteristics. 

Figure  1  shows  two  plots.  In  one  (the  circles),  the  size  dependence  of  the  total 
decomposition  probability  is  given.  In  this  plot,  the  ratio  of  the  sum  of  the  mass 
peak  intensities  of  all  the  daughter  ions  formed  from  evaporation  from  a  certain 
cluster  to  the  sum  of  the  mass  peak  intensities  of  all  the  daughter  ions  and  the 
remaining  parent  cluster  ion  is  given  as  a  function  of  the  cluster  size.  This  is  taken 
as  a  measure  of  the  size  dependence  of  the  total  evaporation  probability.  In  the 
other  plot  (the  triangles),  the  size  dependence  of  the  total  cluster  yield  (i.e.,  the 
relative  intensity  of  the  mass  peaks  of  the  different  parent  clusters  in  the  mass 
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spectrum)  is  given.  It  is  significant  to  observe  the  "anticorrelation"  between  the  two 
plots.  Clusters  that  have  high  yield  in  the  mass  spectrum  (e.g.,  the  ones  with  magic 
numbers  6,9,13)  also  show  low  evaporation  probabilities.  This  strongly  suggests  that 
clusters  that  are  stable  towards  evaporation  indeed  show  relatively  high  mass  peak 
intensities. 

Figure  2  shows  the  size  dependence  of  the  formation  probability  of  the 
daughter  cluster  ions  [Cs(CsI)m]'‘^  by  the  loss  of  either  one  (circles)  or  two  (given  by 
x)  Csl  molecules  from  a  parent  ion  cluster  in  which  n=m  + 1  and  m+2  respectively. 
Again,  the  relative  mass  peak  intensity  of  a  daughter  ion  [Cs(CsI)jjj]‘‘^  to  the  sum  of 
the  mass  peak  intensities  of  a  parent  ion  and  all  the  daughter  ions  is  taken  as  a 
measure  of  the  formation  probability  of  the  daughter  ion  cluster.  This  figure  again 
shows  that  whenever  m  is  a  magic  number  (e.g.  6,9,13...)  the  probability  of  formation 
of  this  daughter  ion  by  evaporation  is  relatively  high. 

Both  of  the  above  two  types  of  results  strongly  support  the  proposal  that  the 
relatively  high  mass  peak  intensities  of  clusters  with  magic  numbers  result  from  their 
low  evaporation  probabilities  as  well  as  their  high  formation  probabilities  from  larger 
clusters  by  evaporating  one  or  more  Csl  molecules. 

C.  Is  the  Evaporation  Process  Statistical?^ 

Does  the  evaporation  of  one  or  two  Csl  molecules  occur  statistically  or  not? 
Does  it  occur  from  an  energy  equilibrated  cluster?  In  our  studies  on  evaporation 
probability,  we  have  observed^^  a  good  anticorrelation  between  the  calculated  beat 
of  evaporation  and  the  evaporation  probability.  This  suggested  statistical  evaporation 
processes  with  no  reverse  activation  energies.  This  conclusion  is  not  surprising  as  the 
condensation  process  (i.e.,  the  reverse  reaction)  represents  an  ion-molecule  reaction, 
known  to  have  zero  activation  energy  in  general. 

Klotz,  in  a  number  of  important  papers,^  discussed  the  evaporation  process 
from  clusters.  The  results  of  his  ev^orative  ensemble  model^  (EEM)  for  statistical 
evaporation  can  be  used  to  obtain^  the  following  relation: 

e^/A  =3/(2Y^)-3kB/(4c(2n+l)]  (1) 

where  is  the  average  kinetic  energy  release;  A  is  the  heat  of  evaporation;  kg  is  the 
Boltzmann  constant;  y  jj  =  Gspan  parameter  for  a  given  microcanonical 

rate  constant  k  (can  be  taken^  as  22  for  k  ~  lO^sec'^,  the  rate  constant  estimated  in 
the  first  field-free  region  in  our  experiment);  and  c  is  the  heat  capacity  for  each 
monomer  unit  (i.e.,  for  each  atom  or  ion).  T  is  the  arithematic  mean  of  the 
temperature  of  the  cluster  at  the  equilibrium  geometry  T  and  that  at  the  transition 
state  T*. 

Temperature  is  used  in  the  EEM  as  an  index  to  the  internal  energy  E  of  a 
cluster  that  evaporates  in  a  given  time  window  (or  with  a  specified  microcanonical 
rate  constant  k),  as  defined  by  E  =  CT  and  T’  =  T  -  A  /C  where  C  is  the  heat 
capacity  of  the  cluster  [i.e.,  C  =  c(2n  +  1)].  If  the  Gspan  parameter,  is  constant 
and  independent  of  the  cluster  size,  Cj/A  goes  as  l/2n+ 1  for  statistical  evaporation. 
y,j  is  proportional  to  A/T  which  has  the  functional  form  of  AS^^p^.  Thus  its 
constancy  for  the  evaporation  of  molecules  from  clusters  resembles  somewhat 
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Trouton's  Rule  for  the  evaporation  of  liquids  or  the  sublimation  of  solids. 

The  general  form  of  equation  (1)  above  can  be  derived  in  a  simple  manner 
as  follows.  The  average  kinetic  energy  released  C2  can  be  expressed  as  follows: 


(2) 

where  E,  A  and  N  are  the  internal  energy  of  the  cluster;  its  heat  of  evaporation  and 
its  internal  degrees  of  freedom,  respectively.  For  a  [CsCCsI)^]''^  cluster,  E=cT(2n+ 1) 
(where  C  is  a  constant  and  T  is  the  absolute  temperature)  and  N  =  3  (2n  +  1)  -  6 
"  3(2n+ 1).  Substituting  for  E  and  N  in  equation  (2)  gives: 

e  =  —  A 

3  '  3(2n+1) 

(3) 

A  3a  3(2n^1) 

(4) 

if  —  is  constant  and  size  independent,  then  (A)  becomes 


€  1 

—  =  constant - 

A  3 


1 

(2m1) 


(5) 


which  is  of  the  same  form  as  equation  (1). 

Figure  3  shows  the  high  resolution  mass  peak  spectra^  of  the  parent  cluster 
ion  [Cs(CsI)g]''’  (bottom),  and  two  of  its  daughter  ions  formed  by  the  evaporation  of 
one  and  two  Csl  molecules  in  the  field  free  region.  For  the  following  decomposition: 
AB^  -  +  B 

the  average  kinetic  energy  release  can  be  calculated  for  Gaussian  peak  shape  (as 
observed  in  our  spectra)  from  the  following  equation:^^ 


c,= 


16/77 


B 


\2 


/  A  \/ 


(6) 

where  is  the  acceleration  voltage  at  the  peak  center  (zero  kinetic  energy  release) 
which  is  =  8.1  keV  in  our  experiment;  m  represents  mass,  A  and  are  the 
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width  at  0.223  of  the  maximum  for  each  of  the  daughter  and  parent  mass  peaks 
measured  in  eV,  respectively.  The  subtraction  in  brackets  corrects  for  the 
instrumental  broadem'ng  as  well  as  the  velocity  and  the  angular  distribution 
broadening  due  to  the  parent  ion  clusters  prior  to  evaporation. 

The  heat  of  evaporation,  A ,  can  be  calculated  by  use  of  the  theoretical  values 
of  the  formation  energies  of  the  clusters^  and  that  of  the  Csl  monomer  or  the 
dimer^.  If  one  uses  these  calculated  values  of  A  and  the  values  of  e,  determined 
from  the  observed  high  resolution  mass  peaks  and  equation  (6)  for  different  clusters, 
one  can  calculate  e^/A  values  for  the  evaporation  of  a  monomer  or  a  dimer  (since 
our  results  discussed  below  suggest  that  the  loss  of  two  Csl  molecules  results  from 
a  fission  of  a  Csl  dimer).  The  plot  of  the  e^/A  values  for  different  clusters  is  given 
in  Figure  4  for  the  evaporation  of  a  Csl  monomer  (circles)  and  the  dimer  (crosses). 
The  theoretical  curve  obtained  by  using  equation  (1)  for  statistical  evaporation 
predicts  the  solid  curve  for  =  22  and  c  =  3kB.  The  agreement  is  indeed  very 
good  and  strongly  suggests  a  statistical  evaporation  of  either  a  monomer  or  the  dimer 
from  these  ionic  clusters.  The  points,  outside  the  predicted  lines  except  for  n  =  3 
-»  2,  can  be  reproduced  (see  Figure  5)  if  we  assume  that  the  temperature  T  instead 
of  Yjj  is  constant.^  This  seems  to  indicate  that  clusters  of  different  sizes  under 
evaporation  in  a  given  time  window  possess  the  same, amount  of  energy  per  each 
monomer  unit.  Although  the  assumption  of  a  constant  T  is  in  a  slight  deviation  from 
the  EEM,  the  much  better  agreement  obtained  with  the  experiments  shown  in  Figure 
5  supports  the  validity  of  this  assumption  and  the  similar  dynamical  characteristic 
"statistical"  nature  of  the  evaporating  clusters  of  differrent  sizes  (except  for  n  =  3  - 
2).  The  large  deviation  observed  for  n  »  3  (see  Figure  5)  might  suggest  either  a 
nonstatistical  evaporation  or  else  an  evaporation  with  large  reverse  activation  energy 
(i.e.,  the  actual  evaporation  process  has  an  activation  energy  larger  than  the 
thermodynamic  heat  of  evaporation).  The  latter  is  the  most  probable  cause  of  the 
deviation  and  results  from  large  energy  of  nuclear  reorganization  in  the  evaporation- 
condensation  process  involving  this  cluster. 

D.  THE  EVAPORATION  OF  TWO  Csl  MOLECULES:  SEQUENTIAL  LOSS  OF 
MONOMERS  OR  FISSION  OF  A  DIMER?^ 

The  evaporation  from  the  different  selected  cluster  ions  usually  gives  daughter 
mass  peaks  resulting  from  the  loss  of  one  or  two  Csl  molecules  that  are  comparable 
in  their  mass  peak  intensities.  If  the  evaporation  was  sequential,  one  would  expect 
the  evaporation  of  the  second  Csl  molecule  to  be  slower  than  that  of  the  first  one 
and  thus  during  the  time-of-flight  less  formation  of  "granddaughter"  ions.  This  is 
because  the  evaporation  of  the  first  molecule  is  expected  to  cool  off  the  resulting 
daughter  ion  cluster  from  which  the  evaporation  of  the  second  Csl  molecule  will  take 
place.  Thus,  as  long  as  the  evaporation  is  statistical  with  no  reverse  activation 
energy,  and  in  the  absence  of  collisions  during  the  flight  time,  one  expects  the  mass 
peak  intensity  of  the  second  daughter  to  be  much  weaker  than  that  for  the  first 
(perhaps,  except  for  the  few  cases  when  the  evaporation  of  the  second  Csl  molecule 
leads  to  a  daughter  ion  with  a  magic  number).  Since  this  is  not  observed,  one  could 
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conclude  that  the  evaporation  of  two  Csl  molecules  occurs  via  a  fission  of  a  (Csl)2 
molecule.  Below  we  present  even  more  quantitative  results  that  substantiate  this 
conclusion. 

a.  Fmm  the  Values  of  the  Kinetic  Ener^  Released  in  Evaporation: 

In  the  absence  of  collisions,  it  is  expected  that  the  evaporation  of  the  first  Csl 
molecule  would  cool  off  the  first  daughter  ion  formed.  The  evaporation  of  a  second 
Csl  molecule  from  this  cold  daughter  is  thus  expected  to  release  less  kinetic  energy 
than  that  released  during  the  loss  of  the  first  Csl  molecule  from  the  parent  ion, 
except  perhaps  in  some  cases  when  the  daughter  resulting  from  the  second 
evaporation  has  a  magic  number. 

We  have  determined^  the  average  kinetic  energy  release  when  a  mass 
selected  [Cs(CsI)„]''’  parent  ion  evaporates  to  give  [Cs(CsI)q.i]'^  and  [Cs(CsI)j,.2]'''. 
Using  equation  (6),  the  average  kinetic  energy  release  for  the  evaporation  of  one  Csl 
molecule  is  calculated  by  using  quantities  for  AB'*’,  A"^,  and  B  corresponding  to 
[Cs(CsI)J'^,  [Cs(CsI)jj.j]‘^  and  Csl,  respectively.  In  order  to  calculate  the  average 
kinetic  energy  release  for  the  evaporation  of  a  second  Csl  from  the  daughter  formed, 
we  use  in  equation  (2)  the  observed  quantities  for  [Cs(CsI)^.j]''',  [Cs(CsI)j,.2]^  and 
Csl  for  AB"^,  A"^  and  B,  respectively.  However,  if  we  want  to  calculate  the  average 
kinetic  energy  release  in  the  evaporation  of  a  dimer  directly  from  the  parent,  we  use 
in  equation  (2)  the  quantities  for  [Cs(CsI)J'*’,  (Cs(CsI)n.2]‘^  and  (Csl)2  for  AB'^,  B"^ 
and  A,  respectively.  The  results^  of  these  three  different  calculations,  i.e.,  the 
evaporation  of  the  first  Csl,  and  the  evaporation  of  the  two  Csl  assuming  sequential 
evaporation  and  assuming  the  fission  of  a  dimer,  strongly  indicates  that  the  dimer 
fission  is  the  most  probable.  The  sequential  evaporation  mechanism  gave  average 
kinetic  energy  release  for  the  evaporation  of  the  Csl  molecule  that  is  greater  than 
(or  comparable  to  in  a  few  cases)  that  for  the  first  one  and  in  one  case  gave  a 
negative  number  for  the  value  of  the  kinetic  energy  release  of  the  second  molecule 
evaporation.  Fission  evaporation  gave  values  that  are  consistent  with  simple  physical 
intuition.  These  results  strongly  support  that  the  loss  of  two  Csl  molecules  from  a 
parent  ion  involves  the  fission  of  a  dimer  rather  than  the  sequential  evaporation  of 
two  Csl  monomers. 

b.  From  the  Distribution  of  Kinetic  Energy  Released 

The  distribution  of  the  kinetic  energy  release  in  the  statistical  dissociation 
(evaporation)  process  is  modulated  by  the  long  range  coupling  between  the  departing 
fragments.  ITius  the  shape  of  the  high  resolution  mass  peak  is  expected  to  reflect  the 
type  of  the  long  range  interaction  between,  and  thus  the  structures  of,  the  fragments 
produced. 

In  sequential  evaporation  from  an  ionic  species,  the  long  range  interartion  is 
expected  to  be  an  ion-dipole  coupling  between  the  remaining  ionic  fragment  and  the 
Csl  dipolar  molecule.  If  this  were  the  mechanism  of  evaporating  two  Csl  molecules, 
one  would  expect  the  mass  peak  shapes  of  the  "granddaughter”  ion  [Cs(CsI)jj,2]^  to 
be  similar  to  that  of  the  daughter  ion  [Cs(CsI)j,.j]''’.  In  both  cases,  ion-dipole 
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interaction  will  be  the  dominant  long  range  coupling.  On  the  other  hand,  if  dimer 
fission  dominates  the  mechanism  of  evaporating  two  Csl  molecules,  the  shape  of  the 
mass  peak  of  the  granddaughter  ion,  [Cs(CsI)^.2]^.  wiW  be  very  different  from  that 
of  the  daughter  ion  ([CsfCsI)^.^]''’).  This  is  because  ion-quadrupole  coupling  will  be 
the  long  range  interaction  involved  in  the  process: 

[Cs(Csi)X  -  [Cs(ai)„.2]^  +  (ai)2 
but  ion-dipole  coupling  will  be  involved  in 
[Cs(CsI)J^  -  ICs(CsI)„.il^  +  Csl. 

The  ratio  Wq  22  7^0.50  ,where  Wq22  and  Wq  jq  are  the  mass  peak  width  at 
0.22  and  0.50  of  the  maxima  (after  correcting  for  the  parent  ion  mass  peak  width), 
is  known  as  the  shape  factor  and  is  a  good  indication  of  the  shape  of  the  high 
resolution  mass  peaks,  and  thus  of  the  distribution  of  the  kinetic  energy  released  ir 
the  evaporation  process.  Table  I  lists^  this  value  for  [Cs(CsI)jj.j] and  [Cs(CsI)„.2]^ 
resulting  from  the  evaporation  of  one  and  two  Csl  molecules  from  clusters  of 
different  sizes  (different  values  of  n).  Independent  of  the  cluster  size,  it  is  clear  that 
all  the  mass  peaks  of  the  daughter  ion  resulting  from  the  evaporation  of  one  Csl 
molecule  have  similar  shapes,  with  a  shape  factor  of  2.30  ±  0.05.  Similarly,  all  mass 
peaks  of  daughter  ions  resulting  from  the  evaporation  of  two  Csl  molecules  have 
similar  shapes,  with  a  shape  factor  of  2.16  ±  0.03,  again  independent  of  n.  This 
result  strongly  confirms  that  the  dissociating  partner  of  the  final  daughter  ion 
resulting  from  the  loss  of  two  Csl  molecules  ([C^(CsI)n.2]''^)  does  not  interact  with 
a  Csl  molecule  (which  is  the  case  in  the  formation  of  the  first  daughter  ion 
([Cs(CsI)j,.2]‘'’).  Thus  it  must  involve  the  interaction  with  a  dimer,  which  has  a  zero 
dipole  moment.  The  values  of  the  shape  factor  indeed  suggest  that  the  loss  of  two 
Csl  molecules  involves  a  much  shorter  range  interaction  than  the  loss  of  one  Csl 
molecule,  which  is  consistent  with  an  ion-quadrupole  and  not  an  ion-dipole  type 
interaction. 

The  question  is  why  would  the  process  of  fission  of  a  dimer  dominate  the 
mechanism  of  evaporating  two  Csl  molecules?  Sequential  mechanisms  have  indeed 
been  observed  to  dominate  the  evaporation  from  some  Cu  and  A1  metal  cluster 
ion.  ’  There  are  several  reasons  for  this  observation:  1)  The  loss  of  one  Csl 
molecule  reduces  the  temperature  of  the  resulting  daughter  ion  cluster  sufficiently 
to  reduce  the  rate  of  evaporating  the  second  Csl  molecule  so  that  it  is  not  observed 
in  the  time  window  remaining  for  it  during  its  flight  to  the  detector.  2)  The  Csl 
dimer  is  a  stable  species,  which  makes  the  heat  of  evaporation  (and  thus  the 
evaporation  activation  energy)  for  the  dimer  much  smaller  than  that  for  the  sum  of 
two  sequential  evaporations  of  two  monomers.  It  is  actually  not  different  from  the 
heat  of  evaporation  of  one  monomer.  This  makes  the  evaporation  of  the  dimer 
compete  well  with  the  evaporation  of  the  monomer  from  the  parent  cluster  ion.  This 
is  the  reason  for  the  observed  comparable  intensity  of  the  mass  peaks  for  [Cs(CsI)^. 
j]'*’  and  [Cs(CsI)n.2]'*^  that  results  from  evaporation  of  one  and  two  Csl  molecules 
from  the  parent  ion  cluster  [Cs(CsI)j,]'^  (except  for  cases  when  clusters  of  magic 
numbers  are  produced). 
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Table  I.  The  Band  Shapes  (as  measured  by  the  ratio  of 

(WidthiQ  22)V(Widthg^)^)  of  the  Corrected  Metastable 
Peaks  of  [(Cs(CsI)j,]‘'‘  after  Loss  of  One  or  Two  Csl 
Molecules  by  Evaporation. 


n 

monomer 

dimer 

2 

2.25 

- 

3 

2.26 

2.14 

4 

2.31 

- 

5 

2.40 

2.17 

6 

2.30 

2.16 

7 

2.30 

2.19 

8 

2.28 

2.19 

9 

2.29 

2.18 

10 

2.24 

2.11 

11 

2.31 

2.16 

12 

2.32 

2.21 

13 

2.26 

2.16 

2.30  ±  0.05 


average 


2.16  ±  0.03 
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Figure  Captions 

Figure  1:  Anticoirelation  between  the  obserwd  intensity  of  the  different  mass  peaks  in  the  mass 

spectrum  of  the  sputtered  [Cs(CsI)„]'*^  dusters  (a)  and  the  evaporation  probability  of 
the  mass  selected  dusters  (•)  The  latter  is  determined  by  allowing  mass  selected 
duster  ions  to  evaporate  dur^  their  flight  to  the  detetSor  and  then  calculate  the 
relatiwe  mass  peak  intenaty  of  all  the  daughter  ions  formed  to  the  total  intensity  of 
these  mass  peaks  and  that  of  their  remaining  parent  duster.  This  anticorrelation 
suggests  that  magic  numbers  observed  in  the  mass  spectrum  result  from  the  stability 
of  the  corresponding  dusters  toward  evaporation. 

Figure  2:  Size  dependence  of  the  probability  of  formation  of  daughter  ions  from  sputtered  mass 

selected  parent  [Cs(CsI)J'^  dusters  by  evaporation  of  one  (solid  lines)  and  two 
(broken  lines)Csl  molecules.  The  fact  that  daughter  ions  with  magic  numbers  (m  = 
6,  9,  13)  haw  the  highest  formation  probability  suggests  that  dusters  with  magic 
numbers  in  the  mass  spectrum  have  relatively  higher  intensity  because  they  don't 
evaporate  to  make  others  (as  shown  in  Figure  1)  and  also  because  other  dusters  of 
higher  mass  evaporate  to  make  them. 

Figure  3:  IVpical  results  of  the  high  resolution  mass  peaks  obtained  from  the  evaporation  of  a 

mass  selected  parent  duster  ion  (e.g.,  [Cs(CsI)5]'^)  to  give  the  first  and  second 
daughter  ions  ([Cs(CsI)5]''’ and  [Cs(CsI)4]‘*^,  respectively).  From  the  width  of  these 
mass  peaks  at  0.22  of  the  maximum  intensity,  the  average  kinetic  energy  released  upon 
the  evaporation  of  one  or  two  (by  sequential  or  fission  type  processes)  Csl  molecules 
can  be  determined  (see  text).  The  band  shape  (see  Table  I)  of  [Cs(CsI)5)*  and 
[Cs(CsI)4]'^,  which  reflects  the  kinetic  energy  distribution  and  thus  the  mechanism  of 
evaporation  of  one  and  two  Csl  molecules,  is  diflerent.  This  strongly  suggests  that  the 
evaporation  of  two  Csl  ion  molecules  is  not  sequential  and  that  the  kinetic  energy 
distribution  of  [Cs(CsI)j]''^  is  modulated  by  ion-dipole  interaction  while  that  of 
(Cs(CsI)4]''’  is  modulated  by  ion-quadrupole  (due  to  the  neutral  dimer)  interaction. 

Comparison  of  experimental  data  and  predicted  (given  by  solid  curve  and  assuming 
statistical  evaporation)  ratio  of  the  average  kinetic  energy  release  to  the  theoretical 
heat  of  evaporation  for  a  Csl  monomer  (o)  and  dimer  (x).  The  predicted  curve  is 
obtained  by  using  Eq.  (1)  with  y^.  =  ZZ  and  c  -  3kB.  The  agreement  suggests 
statistical  evaporation.  The  points  outside  the  predicted  line  can  be  explained  (see 
text). 

Comparison  of  the  measured  (0)  and  the  predicted  {••)  < ,  /  A .  The  predicted  curve 
is  obtained  by  using  Eq.  (1)  assuming  T  is  constant  (T  =  790K)  and  c  «  3kg. 
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